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Frame-independence of the Inhomogeneous Mixmaster Chaos via Misner-Chitre-like 

variables 
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We outline the covariant nature,with respect to the choice of a reference frame, of the chaos 
characterizing the generic cosmological solution near the initial singularity, i.e. the so-called inho- 
mogeneous Mixmaster model. Our analysis is based on a "gauge" independent ADM- reduction of the 
dynamics to the physical degrees of freedom. The out coming picture shows how the inhomogeneous 
Mixmaster model is isomorphic point by point in space to a billiard on a Lobachevsky plane. Indeed, 
the existence of an asymptotic (energy-like) constant of the motion allows to construct the Jacobi 
metric associated to the geodesic flow and to calculate a non-zero Lyapunov exponent in each space 
point. The chaos covariance emerges from the independence of our scheme with respect to the form 
of the lapse function and the shift vector; the origin of this result relies on the dynamical decoupling 
of the space-points which takes place near the singularity, due to the asymptotic approach of the 
potential term to infinite walls. At the ground of the obtained dynamical scheme is the choice of 
Misner-Chitrlike variables which allows to fix the billiard potential walls. 

PACS numbers: 04.20.Jb, 98.80.Bp-Cq 



BASIC STATEMENTS 

The homogeneous and isotropic Friedmann-Lemaitre- 
Robertson- Walker (FLRW) metric provides a valuable 
framework to describe the history of the Universe up 
to its very early stages of evolution. Indeed the very 
good agreement between the light elements abundances 
predicted for the primordial nucleosynthesis and the ob- 
served one allows to extrapolate backward in time the 
FLRW dynamics up to 10~ 2 — 10~ 3 s 0; furthermore, 
recent observations of the Cosmic Microwave Background 
(CMB) radiation HQ, suggest that an inflationary sce- 
nario took place and that the Universe was homogeneous 
and isotropic on the horizon scale up to O(10 _32 )s. In 
spite of such an experimental evidence in favor of the 
FLRW model, there are mainly two well-grounded rea- 
sons to believe that, when the Universe temperature was 
above the Grand Unification scale (O(10 15 GeV)) and 
below the Planck mass (O(10 19 GeV)), the Universe was 
appropriately described by a generic inhomogeneous cos- 
mological model 0,0,0; indeed we have to stress that: 

1. as shown in |7J the Universe in an expanding picture 
is stable with respect to tensorial perturbations; in 
fact the amplitude of a gravitational wave decays 
like the inverse of the scale factor 0, 0] . Therefore 
reversing the expanding behavior into a collapsing 
one, the homogeneity and the isotropy of the Uni- 
verse become unstable with respect to wave like 
perturbations. In particular in @ is shown how 
a Bianchi type IX model far from the singularity 
can be represented (in term of exact solution) by 



a closed FLRW model plus small gravitational rip- 
ples; instead, close to the singularity this model 
has a fully developed anisotropy which results into 
a chaotic behavior. This issue is valuable because 
the Bianchi type IX cosmology presents features 
extensible (point by point in space) to the Generic 
Cosmological Solution. 

2. It is commonly believed 0,E3 that the actual Uni- 
verse came out from a quantum regime which was 
fully developed during the Planckian era and ap- 
proached a classical limit (see [l]}) only in a later 
stage of evolution. Since the early Universe con- 
tained more than a single horizon (see also [T3|L 
then during the Planckian era the metric and topol- 
ogy quantum fluctuations had to take place inde- 
pendently over two causally disconnected regions, 
and hence the survival of any global symmetry is 
prevented. In view of this, we infer that the quan- 
tum behavior of the early Universe has to be prop- 
erly analyzed only in terms of a generic inhomoge- 
neous model. 

With respect to these two points, it is worth noting that 
the bridge between the generic evolution and the FLRW 
dynamics is provided on the horizon scale just by an in- 
flationary scenario, as outlined in [l2T |. 

As well known, the generic cosmological solution Q 
is characterized near the Big-Bang by a chaotic evolu- 
tion which reduces the space-time to the structure of a 
foam 0, 0. When approaching the cosmological singu- 
larity, the space points dynamically decouple and a time 
evolution resembling the so-called Mixmaster dynamics 



of the Bianchi VIII and IX models takes place indepen- 
dently in each of them [l3l 1 1 ol | (for further discussions on 
inhomogeneous cosmological models see [la,ll2|); from a 
physical point of view, space neighborhood is here consid- 
ered at horizon size (for recent discussions on Mixmaster 
covariance see|la llSj |20|). 

In recent years, the chaoticity of the homogeneous Mix- 
master model has been widely studied in the literature 
(see [H,|il,|23,|2ll|23,|23|) in view of understanding the 
features is of its covariant nature. 

Two convincing arguments, appeared in [jol l22l| . support 
the idea that the Mixmaster chaos (described by the in- 
variant measure introduced in |23l l24j) remains valid in 
any system of coordinates. 

The main issue of the present work is to show that the 
property of space-time covariance can be extended to the 
inhomogeneous Mixmaster model. In Section 2 we pro- 
vide a gauge independent analysis (i.e. independent of 
the choice of the lapse function as well as of the shift 
vector) for the dynamics of the gravitational degrees of 
freedom. In Section 3 we discuss the asymptotic behavior 
of the potential term associated with the Ricci 3-scalar, 
showing how it can be modeled in terms of a potential 
wall surrounding a (dynamically) closed domain. In Sec- 
tion 4 in view of the existence of an energy-like constant 
of motion, we construct the Jacobi metric associated with 
the resulting billiard on a two dimensional Lobachevsky 
plane point by point in space, and finally we calculate the 
covariant non-zero Lyapunov exponent. The concluding 
remarks presented in Section 5 are devoted to outline 
how the approximation used for the potential term is 
completely self-consistent being dynamically induced. It 
is worth noting how we show that for any choice of lapse 
function and the shift-vector, taking Misner-Chitre like 
variables, the system results to be chaotic; however, in 
this sense, our analysis states the independence of the 
chaos on the choice of the space-time coordinates, but 
the questions of covariance remains open when using dif- 
ferent configuration coordinates because the Lyapunov 
exponent is not a good indicator for all of them. 



HAMILTONIAN FORMULATION 

A generic cosmological solution is represented by a 
gravitational field having available all its degrees of free- 
dom and, therefore, allowing to specify a generic Cauchy 
problem. 

In the Arnowitt-Deser-Misner (ADM) formalism, the 
metric tensor corresponding to such a generic model takes 
the form 



dT 2 



N 2 dt 2 - la p{dx a + N a dt){dx p + N p dt) (1) 



where N and N a denote (respectively) the lapse function 
and the shift-vector, ~f a p (a, j3 = 1, 2, 3) the 3-metric ten- 
sor of the spatial hyper-surfaces E 3 for which t — const, 



being 

7a/3 = e^5 ad O a b O d c d a y b dpy c , a, b, c, d,a,/3= 1, 2, 3, 

(2) 

while q a = q a (x, t) and y b = y b (x, t) six scalar functions, 
and Ol = 0%(x) a SO(3) matrix. By the metric tensor 
, the action for the gravitational field is 

S = [ dtd 3 x( Pa d t q a +Il d d ty d -NH-N a H a ) 



where 



(3) 



V' a b 



H a = U c d ay c + Pa d a q a + 2p a {0- l ) b a d a O a b ; (5) 

in j3j and JHJ) p a and Tld are the conjugate momenta 
of the variable q a and y b respectively, and the is 
the Ricci 3-scalar which plays the role of a potential term. 

The ten independent components of a generic metric 
tensor are represented by the three scale factors q a , the 
three degrees of freedom y a , the lapse function N and 
the three components of the shift-vector N a ; it is worth 
noting that, by the variation of the variables p a , Tl a in 
the action ©, the relations: 



d t y d = N a d a y d 



N = 



(6) 



(7) 



take place. 



ADM- REDUCTION OF THE DYNAMICS 

We use the Hamiltonian constraints H = H a = for 
the reduction of the dynamics to the physical degrees 
of freedom; from Q and JSJ, we note that the super- 
momentum constraints can be diagonalized and explicitly 
solved by choosing the function y a as special coordinates, 
i.e. taking the transformation rj = t, and y a — y a (t,x); 
in fact starting by JU and © we get the expression 
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dq a 



-i, c dO a c 



(8) 



Furthermore, in the new coordinates we have 

q a (t, x)^q a {n,y) 
p a (t,x) -^p' a (r),y) =p a (r),y)/\J\ 

JL £>2L._£L 4_ JL 

dt dt 0y b ^ dr) 

_d ¥ dy b d 

v dx a dx a dy b 



(9) 



where |J| denotes the Jacobian of the transformation. 
The first relation holds in general for all the scalar quan- 
tities, while the second one for all the scalar densities; 
hence the action Q rewrites as 



S= f d V d 3 y ( Pa d v q a + 2 Pa (0- 1 ) c a d v O a c - NH) 

JSC 3 ) xSR 



(10) 



THE POTENTIAL WALL AND THE REDUCED 
VARIATIONAL PRINCIPLE 

The potential term appearing in the super- 
Hamiltonian reads, in obvious notation, as: 

U = tjL (3) i? - £ XlD 2 ^+J2 D^O (dq, (dqf, y, , 



b^c 



where 



# = exp£<f, 



Qa 



kj 
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(11) 

(12) 

(13) 
(14) 



Assuming the functions y a {t, x) smooth enough (which 
implies by ©,0 that the coordinates system is smooth 
"itself), then all the gradients appearing in the potential 
U are regular. Indeed this notion of regularity is not to 
be intended in absolute sense; in fact what really mat- 
ter here is not that the gradient increase but simply that 
their behavior is not so strongly divergent to destroy the 
billiard representation (see next paragraph). In [5j it was 
shown that the spatial gradients increase logarithmically 
in the proper time along the billiard's geodesic and there- 
fore result to be of higher order. Thus, as D — > [3(| 
the spatial curvature diverges and the cosmological 
singularity appears; in this limit, the first term of the 
potential U dominates all the remaining ones and can be 
modeled by the potential wall 



being 



e(x) 



+oo if x > 0, 
if a; < 0. 



(15) 



(16) 



By ()15f) the Universe dynamics evolves independently 
in each space point; the point-Universe can move only 




FIG. 1: r Q (£,9) 



within a dynamically-closed domain Tq (see figure l|T]l) 
(indeed the three corners at infinity are open but they 
correspond to a set of zero measure in the space of ini- 
tial conditions). Since in Tq the potential U asymptoti- 
cally vanishes, near the singularity we have dp a /di] = 0; 
and then the term 2p a (0^ 1 ) c a d, 1 < ^ in (|ll)fl behaves as 
an exact time-derivative and can be ruled out of the 
variational principle. The ADM reduction is completed 
by introducing the so-called Misner-Chitre like variables 
|5j, 1221 |25j, |26j as 



e T W£ 2 - 1 (cos 
= ^W? ~ l(cos0 



V3 sin 8 
- \/3sin( 



(17) 



1 cos ( 



The way in which the anisotropy parameters Q a Ijl^jl 
are rewritten is an important feature of these variables 
since we easily get the r-independent expressions: 



Qi = i ^ 1 (cos 9 + \/3sm( 



Qi = \ ^7 1 ( cos - V3sin( 



(18) 



Qa 



2ve 



When expressed in term of such variables the super- 
hamiltonian constraint can be solved in the domain Tq: 



and the reduced action reads as 

SSr Q = 8 J drjSyip^d,^ + pgd v 9 - ed v r) = . 



(19) 



(20) 



By the asymptotic limit l|15f) and the Hamilton equations 
associated with (|20|) it follows that e is a constant of 
motion, i.e. de/drj — de/d?] = =>■ e = E(y a ). 



THE JACOBI METRIC AND THE LYAPUNOV 
EXPONENT 

Being e a constant of motion, the term ed v r = 
E(y a )d v T in H2()(l behaves as an exact time derivative; 
hence the variational principle rewrites: 



d 3 y( Pi d^+ P gd9) =0, 



(21) 



coupled with the constraint l|19(l . This dynamical scheme 
allows to construct the Hamilton-Jacobi metric [2?} cor- 
responding to the dynamical flow. Indeed for each point 
of the space it can be reproduced the same analysis de- 
veloped in [2^] for the homogeneous Mixmaster model; 
in particular all the spatial gradients are dumped and 
the space points dynamically decouple in the asymptotic 
limit to the singularity. In fact, in each space point, the 
system dynamics is replaced by a geodesic flow 6 J ds = 
0, with 



ds 2 = E\y a ) 



(e 1) 



+ (e 2 - !)<& 



(22) 



corresponding to the Jacobi line element. The Jacobi 
metric is valid independently in each of such point (here 
the space coordinates behaves like external parameter) 
since the evolution is spatially uncorrelated. The Ricci 
scalar takes the value R = —2/E 2 , hence such a metric 
describes a two-dimensional Lobachevsky space; the role 
of the potential wall (|15(l consists of cutting a closed 
domain Tq on such a negative curved surface. Thus, 
summarizing, the system obtained is isomorphic to a 
billiard on a Lobachevsky plane. 

A precise information about the dynamical stability of 
the geodesic flow associated with the line element <|22[) 
arises by the calculation of the Lyapunov exponent. Let 
us project the connecting vector between two close 
geodesies on a Fermi basis w^}, (fi — 1,2); the 
geodesic vector u M is taken in the form 



d£ d9\ 
ds' ds J 



ve 



E 



1 sin0(s) 
— cos mis) 



h E^e 



1 



(23) 

where s denotes the curvilinear coordinate, while 4>{s) 
is an angular variable (0 < </> < 27r), whose dynamics 
is obtained by requiring that the geodesic equation is 
verified, i.e. 



#0) 

ds 



Eve 



1 



sin (f>(s) . 



(24) 



The vector is determined by the property of the Fermi 
basis to be orthonormal, and it reads explicitly 



ur 



Vt 2 - 1 • j cos< 
■ sin & 



E 



Eve^r 



(25) 



Projecting the connecting vector Z^ over the Fermi basis 
defined above 



Z» = Z u {s)u» + Z w {s)w\ 



(26) 



from the geodesic deviation equation we get that the dy- 
namics is described by the following equations 



(27) 




ds 2 E 2 

The solution for the system (|27|l is given by 




(28) 



The value of E given by the constraint (|19fl and involved 
in the line element l|22|l is determined by the initial condi- 
tions and cannot vanishes. By the first of solutions l|28|) 
no geodesic deviation takes place along the geodesic vec- 
tor (as expected); instead from the second solution we 
get a non-zero Lyapunov exponent of the form 



X(y a ) = limsup 



\n(Zl + (dZ w /ds) 2 ) = 1 

2s E(y») 



> 0. 



(29) 

For the validity of this analysis we have to verify that in 
the limit to the initial singularity the curvilinear coordi- 
nate s approaches infinity; repeating the same procedure 
of [23, formulas (31) appearing there, in our inhomoge- 
neous case is replaced by 



ds 
drj 



E{y a )d n T 



(30) 



i.e. (being f(y a ) a generic function of the space coordi- 
nates) 



E(y a )r + f(y a ) 



lim s 

T— »00 



(31) 



this ensures that the curvilinear coordinate s behaves in 
the appropriate manner in the limit toward the singular- 
ity (t -> 00). 

Hence the chaoticity of the inhomogeneous Mixmaster 
dynamics is ensured by Tq to be a closed domain |3l| : 
the covariance of such a description follows from the in- 
dependence of the Lyapunov exponent with respect to 
the lapse function and the shift vector. In fact, in © 
and 10 N and N a are fixed (in turn) by choosing the 
form of the quantities y a and the latter can be generic 
functions subjected only to the condition to be smooth 



enough. Eq. i|29[l provides the form of the Lyapunov ex- 
ponent in the whole space domain but we stress how its 
value depends on the choice of Misner-Chitre like vari- 
ables; the independence of our scheme on the shift vector 
is ensured by the asymptotic behavior of the potential 
term, but to get e as constant of motion, allowing the 
Jacobi metric representation, we need Misner-Chitre like 
variable. 

The covariance of our picture is equivalent to the co- 
variance of the inhomogeneous Mixmaster chaos because 
it is well known [22I Ejf that the obtained billiard has 
stochastic properties (see also H^). In fact the negative 
curvature of the Lobachevsky plane makes unstable the 
geodesic flow; the potential walls have the role of replac- 
ing a given geodesic with a different one (whose tangent 
vector is related to the previous one by a reflection rule 
|28|) and as we will show in Section 6 their structure will 
influence the chaotic properties of the system dynamics. 

To better characterize the chaoticity of the obtained 
billiard, we show that in each point of the space our sys- 
tem admits an invariant measure which in the present 
variable is uniform over the admissible phase space. 

From the point of view of statistical mechanics, such 
a system admits, point by point in space, an "energy- 
like" constant of motion which corresponds to the kinetic 
part of the ADM Hamiltonian e = E{y a ). The point- 
universe, randomizing within the closed domain Tq, is 
represented by a dynamics which allows for an ensemble 
representation; in view of the existence of the "energy- 
like" constant of motion, the system evolution is appro- 
priately described by microcanonical ensemble. There- 
fore the stochasticity of this system is governed by the 
Liouville invariant measure 



dp oc 8 {E(y a ) - e) d^dOdp^dpe 



(32) 



where S (x) denotes the Dirac functional. 

Since the particular value taken by the function e (e = 
E{y a )) cannot influence the stochastic property of the 
system and must be fixed by the initial conditions, then 
we must integrate (in functional sense) over all admissible 
form of e. To do this it is convenient to introduce the 
natural variables (e,tp) in place of (p£,pe) by 



Pi 



: COS (,5 



= e VC 2 - lsin<£, 



(33) 



where < ip < 2n. By integrating over all functional 
forms of e,we removes the Dirac delta functional, which 
leads in each point of space to the uniform normalized 
invariant measure |32j 



The existence of the above stationary probability dis- 
tribution in Tq, outlines the chaotic properties associated 
to the point-like billiard out coming from our analysis. 

As we stressed, our dynamical scheme relies on the use 
of Misner-Chitre like variables and therefore the covari- 
ance of the Lyapunov exponent is invariant with respect 
to space-time coordinates, but it could be sensitive to the 
choice of configurational variables. Thus the result here 
obtained calls attention to be extended to any choice of 
the configurational variables. In this respect we compare 
our result with the analysis presented in .19] according 
to which, given a dynamical system of the form 



dx/dt = F(x) 



(35) 



d/i(y a ) = d£d0d(p 



8tt 2 ' 



(34) 



then the positiveness of the associated Lyapunov expo- 
nents are invariant under the following diffeomorphism: 
y = </>(x, t),dr = A(x, t)dt, as soon as the four require- 
ments hold: 

1. the system is autonomous 

2. the relevant part of the phase space is bounded 

3. the invariant measure is normalizable 

4. the domain of the time parameter is infinite 

To show that such a covariance criterion is here ful- 
filled, we observe that the variables x can be identified 
with r, £, 6 and 9?, and the time variable with our curvilin- 
ear coordinate s. On the other hand the above diffeomor- 
phism relation in its time independent-form can match a 
phase-space coordinates transformation; then we under- 
line also that: 

1. in the considered asymptotic limit our dynamical 
system is autonomous because its hamiltonian co- 
incides with the constant of motion e = E(y a ) and 
the potential walls are fixed in time. 

2. apart from sets of zero-measure which cannot be 
explored by the system 0, the phase-space Fq is 
a compact domain. 

3. as shown by the above analysis which leads to Ij34(l 
the system admits, in each space point, a normal- 
ized invariant measure over the phase space. 

4. we showed by (|31|l that the curvilinear coordinate 
s admits an infinite domain because the variables 
r G (—00, 00). 

Thus, on the base of ,we can claim that the Lya- 
punov exponent calculated in l|29|l provide an appropriate 
chaos indicator only when the effects of the boundary are 
taken into account in agreement with our discussion of 
the next section. Furthermore such an indicator is covari- 
ant with respect any configurational coordinate transfor- 
mation which preserves the requirements 1-4. Strikingly, 



we have to stress that if we adopt Misner like variables 
|l5| the Lyapunov exponent l|29|> is no longer a good indi- 
cator; in fact the anisotropy parameters (|13|1 in Misner- 
like variables depend not only on but also on a, 

and therefore after the ADM reduction on the curvilinear 
coordinate s. As a consequence, the conditions 1 and 3 
are no longer full filled for this choice because the po- 
tential walls move with time. However for any generic 
transformation of coordinates which involves only £ and 
8, the chaoticity of the Mixmaster model is preserved. 
We conclude this section by observing that if in different 
configurational coordinates the inhomogeneous Mixmas- 
ter would not appear as chaotic, then the stochasticity 
has to be transferred to the coordinate transformation 
which link them to the Misner-Chitre like variable. 



THE ROLE OF THE POTENTIAL WALLS 
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FIG. 2: T Q {u,v) 



Up to now we involve the role of the potential walls in 
order to cut a billiard on the Lobachevsky plane. Here 
we discuss a notion of Lyapunov exponents which include 
either the feature of the geodesic flow cither the structure 
of the bounding potential walls, and we arrive to show 
that even in this more complete framework our system 
is a chaotic one. To this end we are showing below that 
our system meets all the hypothesis at the ground of 
the Wojtkowski Theorem, see [2^. Let us consider the 
following new choice of coordinates on the 2-surface: 



P = (2/1,2/2) 



(1+0 



(cos 9, sin ( 



(36) 



On the basis of i|36|) the new line element and the 
anisotropy parameters read as 



ds 2 = 



AE 2 df 
(l-V 2 ) 2 ' 



y < l; 



[{f+A a f + \-{A a f 



1,2,3 



(37) 



(38) 



being Ai(-V3,V3,0) and ^(1,1,-2) 

Since here we introduce the Poincare model of the 
Lobachevsky plane in the form of the upper-half plane; 
this is reached by using the well-known Poincare variables 



V + b 

2 y ^ -b 

(y + b) 2 



(39) 



where b denotes a point on the absolute (b 2 = 1). In 
terms of these new variables the metric (|36fl takes the 
form 



ds 2 



(dfff 



(40) 



and the available domain \y\ < 1 transforms into the 
half-plane (if ■ b) > while the absolute represents the 
line (77 • b) = 0. Now if we write 



2 

71 











(u-\ 






+ vb 




~2y 







v>0 



(41) 



where b = (0, 1), b = (1, 0), then is easy to verify that in 
term of these coordinates the anisotropy functions have 
the form 



)i(u,v) 



u(u+l)+v 2 

U 2 + U+1 + V 2 



(42) 



This is a very suitable expression for the boundary; 
in fact geodesic on this half-plane are semicircles having 
centers on the absolute (i.e. u(s) = A + R cos s,v(s) — 
Rsms being (A, 0) the coordinate of the semicircle center 
and R the corresponding radius) and rays being perpen- 
dicular to the absolute; the billiard is bounded by the 
geodesic triangle u = 0,u = —1, and (it + 1/2) 2 + v 2 = 
1/4. The new domain is shown in fig 2. 

The billiard has a finite measure, and its open region 
at infinity together with the 2 points on the absolute 
(0,0) and (—1,0) correspond to the three cuspids of the 
potential in figure 1. The tangent field to the geodesic 
flow takes the explicit form: 



f = (-v,u- A) 



(43) 



and the associated matrix of the dynamics M" = d^f is 
constant and orthonormal. 

Referring to the notation of |29j we can easily construct 
in this variables the cocycle of the dynamics and verify 
that theorem 2.2 (pag 151 of this reference) applies; since 
the manifold is connected, has finite volume, and the 
matrix of the cocycle is constant, we have only to find an 



invariant bundle of sectors. This is possible because the 
following two properties hold 

1. because of the constant negative curvature of the 
surface, a one-parameter family of geodesies with 
negative curvature is invariant under the evolution 

2. the bounding potential walls are constituted by two 
straight lines and a semicircle of negative curvature; 
the former ones do not affect the structure of the 
cones during the bounces of the geodesic, while the 
latter one, being a dispersing profile, ensures that 
after reflection against it, the cones will evolve in 
themselves. 

After this discussion we can claim that the largest 
Lyapunov exponent has positive sign almost everywhere. 
The covariance of this result is ensured in view of the dis- 
cussion developed in the previous section based on C[l9|l. 
The analysis of this section complete our analysis about 
the covariant nature of the chaos associated to the bil- 
liard representation of the Mixmaster model. 

CONCLUDING REMARKS 

The main issue of the present work consists of the 
proof that the chaotic behavior singled out by a generic 
inhomogeneous model near a singularity has a covariant 
nature. This result has been obtained by a "gauge" 
independent ADM reduction of the dynamics to the 
physical degrees of freedom, which for the Universe cor- 
respond to the anisotropy degrees, i.e. to the functions 
£ and 6. 

We describe the evolution as independent of the lapse 
function and the shift vector form by adopting the vari- 
ables y a as the new spatial coordinates. However their 
degrees of freedom do not disappear from the problem 
because they are transferred to the SO(3) matrices 
which acquire a dependence on time in the new vari- 
ables; such a dependence is then eliminated from the dy- 
namics by solving the super-momentum constraint and 
using some implications deriving by the approximation 
of y/g^R as an infinite potential wall. The potential be- 
havior l |l 5(1 is crucial for the existence of an " energy-like" 
constant of motion e = E(y a ), and therefore is at the ba- 
sis of the chaos description. In this view, this approxima- 
tion is naturally induced for D — > 0, due to the potential 
structure; the only assumption required is for the func- 
tions y a and 0%(y c ) to be smooth ones, in the sense that 
their presence does not affect the asymptotic behavior 
of the potential term. This restriction is a natural one 
having a good degree of generality, because the smooth- 
ness of the functions y a is ensured by the smoothness of 
the lapse function and the shift vector (see 1 (617(1 1. i.e. 
by the choice of a regular reference frame. The initial 



smoothness of the matrices 0% (when are taken on the 
spatial coordinates x a ) is preserved by the coordinates 
transformation. Concluding, the cosmological meaning of 
this concept corresponds to deal with independent " hori- 
zons" ; the approximation neglecting in 119(1 the potential 
term with respect to the value of " e" is equivalent to re- 
quire the scale of the inhomogeneities to be super-horizon 
sized (see 01111). 
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